Spin-orbit interaction in symmetric wells with two subbands 
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We investigate the spin-orbit (s-o) interaction in two-dimensional electron gases (2DEGs) in quan- 
tum wells with two subbands. From the 8x8 Kane model, we derive a new inter-subband-induced 
s-o term which resembles the functional form of the Rashba s-o - but is non-zero even in symmetric 
structures. This follows from the distinct parity of the confined states (even/odd) which obliterates 
the need for asymmetric potentials. We self-consistently calculate the new s-o coupling strength for 
realistic wells and find it comparable to the usual Rashba constant. Our new s-o term gives rise 
to a non-zero ballistic spin-Hall conductivity, which changes sign as a function of the Fermi energy 
(ef), and can induce an unusual zitterbewegung with cycloidal trajectories without magnetic fields. 

PACS numbers: 72.25.Dc,73.21.Fg, 71.70.Ej, 85.75.-d 



The rapidly developing field of spintronics has gener- 
ated a great deal of interest in s-o coupling in semicon- 
ductor nanostructures For an n-doped zincblende 
semiconductor quantum well with only the lowest sub- 
band occupied, i.e. in a strictly 2D situation, there are 
two main contributions to the interaction of the spin and 
orbital degrees of freedom of electrons. One contribution 
is the Dresselhaus term, which results from the lack of 
inversion symmetry of the underlying zinc-blende lattice 
2j and is to lowest order linear in the crystal momentum 
3|. This linearity is shared by the other contribution 
known as the Rashba term Q , which is due to structural 
inversion asymmetry and can be tuned by an electric gate 
across the well These two contributions can lead to 
an interesting interplay in spintronic systems [||. 

In this Letter we consider yet another type of elec- 
tronic s-o coupling which, as we show, occurs in III-V 
(or II- VI) zinc-blende semiconductor quantum wells with 
more than one subband. We derive a new inter-subband- 
induced s-o interaction which resembles that of the ordi- 
nary Rashba model; however, in contrast to the latter, 
ours is non-zero even in symmetric structures, Fig. 1. 
We self consistently determine the strength of this new 
s-o coupling for realistic single and double wells and find 
it comparable to the Rashba constant, Figs. 2(a)-(b). We 
have investigated the spin Hall effect and the dynamics 
of spin-polarized electrons due to this new s-o term. We 
find (i) a non-zero ballistic spin Hall conductivity which 
changes sign as a function of £p and (ii) an unusual zit- 
terbewegung Q with cycloidal trajectories without mag- 
netic fields, Fig. 3. As derived below, for a symmetric 
well with two subbands our 4x4 electron Hamiltonian is 
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FIG. 1: Square well with its ground-state <fi e (z) and first 
excited-state tp (z) wavefunctions. The new inter-subband- 
induced s-o coupling r\ in Eq. ([2]) is non-zero even in sym- 
metric wells due to the distinct parities of <p e (z) (even) and 
<fio(z) (odd), which yield a non- vanishing matrix element for 
the derivative of the symmetric potential. 



e are quantized energies of the lowest (even) and first 
excited (odd) subbands (corresponding to eigenstates |e) 
and |o)), respectively, measured from the bottom of the 
quantum well, t x,v,z denote the Pauli matrices describ- 
ing the subband (or pseudospin) degree of freedom, and 
a x,y,z are Pauli matrices referring to the electron spin. 
The new inter-subband-induced s-o coupling n is 
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where m* is the effective mass, e± = (s ± s e ) /2, e e and 



where E g and A are the fundamental and split-off band 
aps in the well region Q ; P is the Kane matrix clement 
The parameters S v and <5a denote valence band off- 



2 



sets between the well and the barrier regions [l(|, V(z) 
is the Hartree-type contribution to the electron potential 
and h(z) is the structural quantum- well profile [llj. Note 
that 77 can be varied via external gates, Fig. 2. Next we 
outline the derivation of TL in Eq. (Q]) . 

Kane Hamiltonian. We start from the usual 8x8 Kane 
Hamiltonian describing the s-type conduction and the p- 
type valence bands around the T point 
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where H c is a 2 x 2 diagonal matrix with elements 
p 2 /2mo+V c (r), mo is the bare electron mass, H v is a 6x6 
diagonal matrix with elements p 2 /2itiq+V v (f) — E g for the 
heavy- and light-hole bands and p 2 /2m + Va (t) —E g — A 
for the split-off band, Vi(r) [i = c, v, A] denote arbitrary 
potentials (see below), and H cv = (H vc y is 
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where k = Pk, k = p/H is the electron wave vector, 
k± = k x ± iky, and P = —ih(S\p x \X) / 'mo parameterizes 
the conduction-to- valence band coupling; 15) and \X) are 
the usual periodic Bloch functions at the T point. 

Effective electron Hamiltonian: folding down. The 
Kane Hamiltonian ((3]) acts on an eight-component spinor 
\f f t = ( tp c ip v ) in which the last six components 
represent valence-band states. By eliminating the hole 
components from the Schrodingcr equation 7^8x8 'I' = , 
where e is the eigenenergy, we can fold down this 8x8 
equation into a 2 x 2 effective equation for the conduction- 
band states only: H(e)$ c = [H c +H cv (e — H v )~ H vc ]i[> c , 
ip c is a renormalized conduction-electron spinor. 

s-o in symmetric wells. Applying the above procedure 
to a quantum well, defined by the confining potentials 
lH Vi(r) -► Vi{z) = V(z) + Sih(z), i = c,v, A, we find 



W(e) = H QW + — p_ [ri^ 1 + ^ ,p z ] , 
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where l/m*(z, e) 

(p 2 /2m + V(z) - S v h(z) - E g ) and i] 2 
V(z) — SAh(z) - Eg - A). Equation ^ describes 
an electron in a quantum well (Hqw term) with spin 
orbit interaction (last term) [131 ]. The kinetic-energy 
operators above are complicated due to the position- 
and energy-dependent effective mass m*(z,e). Since E g 
and Eg + A are the largest energy scales in our sys- 
tem, we can simplify §5§ and ((6]) by expanding 1/rji 
and I/772 in the form I/771 = E^fl — [e — p 2 /2mo — 
V(z) + 5 v h(z)]/E g 
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zeroth order rji = 
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8 A h(z)}/ (E g + A) + ...}. To 
Eg, T]2 = Eg + A and H QW = 
V c (z) with (a constant effective mass) 

Since the s-o op- 
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V{E g 
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d z (l/r)i) + d z (l/n 2 ), we need 
to keep the first-order terms in the expansions of 77^ 
and 1J2 1 which yield the leading non-zero contribution to 
We find [ V y+r 1 ~ 1 ,p z ] = [1/E 2 - 

5a/ {Eg + Af]d z h(z). 
Finally, we project this s-o operator into the two lowest 
(spin-degenerate) eigenstates \i) az = \k\\i)\cr z ), (r\k\\i) = 
exp(ifcii -r\\)ipi(z), i = e,o and a z =t,4, of the symmetric 
well (Hqw), Fig. 1. This directly leads to the Ti. in (Q]) 
with the s-o coupling 77 ([2]) [15j |. Note that this new s-o 
interaction is non-zero even in symmetric wells as 77 arises 
from the coupling between the ground state (even) and 
the first excited state (odd), Eq. We can generalize 
H to include the Rashba a and the linearized Dressclhaus 
(3 s-o couplings. Next we determine the magnitude of 77 
(and a, (3) for realistic quantum wells with two subbands. 

Self- consistent calculation of the s-o couplings. We 
consider modulation-doped quantum wells similar to 
those experimentally investigated in Ref. [16]. Our 
wells, however, have two occupied subbands. Similarly 
to Ref. lid], we study cases with constant chemical po- 
tentials [17J. By self-consistently solving Poisson and 
Schrodinger's equations we determine the energy levels 
e e , e and the confined wave functions ipi(z), i = e, o of 
the wells. We then calculate (i) 77 via Eq. ((2]), (ii) a» 
from equations similar to Eq. @ for each subband, and 
(iii) /3i = j c (i\k 2 \i), 7 C is the bulk Dressclhaus s-o param- 
eter [l8| . The structural symmetry of the wells and their 
charge densities can be changed via a gate potential 

Our calculated s-o couplings 77, on, and ft, i — e, o, for 
an InAlAs/InGaAs/InAlAs single quantum well ("sam- 
ple 3" in [l6() are all comparable in magnitude, Fig. 2(a). 
Note that our o; e /ft ratio is consistent with the experi- 
mental one in Ref . (lij . In addition, a e vs. V& here agrees 
well with the experimental data in Fig. 3 ( "tr iang le up") 
of Ref. (see also Fig. 4 in 0). Our a % and ft 
are also consistent with those of Ref. [13] ■ Note that for 
the single well studied here r\ docs not vary appreciably 
with the gate Vb, similarly to ft and as opposed to oti. 

For a double-well structure, on the other hand, we find 
that 77 has a "resonant behavior" , changing by about an 
order of magnitude as Vb is swept across Vb = 0, Fig. 2(b) 
[this may have a dramatic effect on Shubnikov-dc Haas 
measurements] . Vb = corresponds to a fully symmetric 
double well. In contrast to the single-well case, a e and 
a a have opposite signs and undergo abrupt changes in 
magnitudes near Vb = [Fig. 2(b), dashed lines]. Simi- 
larly to the single well case, ft and ft are also essentially 
constant for a double well [Fig. 2(b), dotted lines]. A de- 
tailed account of our results will be presented elsewhere. 
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FIG. 2: Calculated s-o coupling strengths as a function of 
the external gate Vb for realistic wells, (a) For the single 
GalnAs [Irjl . l20l | well studied, the inter-subband coupling n is 
larger than the Dresselhaus j3i and the Rashba at constants 
(i — e,o). Note that |a e | > \at \ and both change sign across 
Vb = (in constrast to /3i and n). (b) For the InSb double 
well considered, n shows a "resonant behavior" about 14 = 
[symmetric configuration, lower- left inset in (b)]. This occurs 
because the subband splitting e — £e reaches a minimum at 
Vb = and the double-well wave functions are very similar 
(though of distinct parities) for Vb ~ 0. This also makes 
a e ~ — a around Vb — 0. Upper-right inset in (b): energy 
dispersions e±(k) [Eq. ([Sj) ] of the symmetric double well. 



Having established that the new s-o coupling -q is sizable, 
in what follows we focus on a fully symmetric well to 
investigate physical effects arising solely from rj. 

Fully symmetric case: eigensolutions. Let us consider 
a two subband well (single or double) described by the 
Hamiltonian 7i in |T]) (we assume a negligible Dresselhaus 
term In the basis {|e)j, |o)j., |o)f, |e)jj H becomes 
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Both the upper-left (U) and lower-right (L) blocks of Ti. 
have eigenvalues 



e±(k) =e k ± hft, 



(8) 



with efc = |^ir + e+, (M7) 2 = ?7 2 fc 2 + e^, and eigenvectors 
|^i>+ = sin (0/2)|e) T + cos (0/2) e^, (9) 

|V>2>£ = COS(e/2)| ) T +8111(^/2)6^16)4, (10) 

|^3)^ = cos^/2)|e) T -sin(^/2)e^| ) i , (11) 
\^) L _ = sin(0/2)|o) T -cos(0/2) e ^|e) i . (12) 



Here. 



{-ky + 



ik x )/k, cos{6) = l/yjl + {vk/e-f, 
and k = A; (sin 4>, — cos 0) (here we drop the "||" in fcy). 
For rjk << 2e_ we can expand £±(fc) in ([5]) and define 
effective masses m* ± = m*/[l ± 2e so /e_], where e so = 
rj 2 m* /2h 2 is the energy scale of the new s-o coupling. 
For the double well of Fig. 2(b), m*_ is reduced by ~ 5% 
compared to the bulk value m* . This could be measured 
via, e.g., cyclotron-resonance experiments [23| . 

Novel Zitterbewegung. The dynamics of electron wave 
packets in wells with s-o interaction exhibit an oscillatory 
motion 0] - the zitterbewegung. For our new s-o interac- 
tion, a wave packet \x) moves according to (xlrff(i)lx) 
where Fjj(t) = U'rU is the position operator in the 
Heisenberg picture [U = exp(—iHt/h)] with components 



XffOt) = 1 <g> lx(0) + 1 <g> 1— t+^tr x ' 

m* n 
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e^r 1 <g> (T y + e_-r 2 



2(M}) 3 

(|) 2 ^ ® (Px<r w +jyr e )] (sin (2^) - 2tl%3) 



and yn(t), obtained from Eq. (|13|) via the replacements 
(p^cr 1 ) 1 ^ (Py,cr y ), o*) >-> (-p^-er 1 ) [i.e. a tt/2 
rotation about the z-axis]. Similar expressions can be 
derived for the spin components o~ l H (t), i = x,y,z 

For simplicity, we evaluate the expectation value of 
th (t) for planes waves ( "wide wave packets" ) . For a spin- 
up electron injected into the lowest subband along the 
y-axis with (group) velocity v g = {KkQ V / rn*)y we find 



Mt)) 
(y H (t)) 



V 2 ko y 
2(M7) 2 



(1 - cos (2^)) , 



ifk 
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2(my 



(sin (2Qt) 



(14) 
2flt) , (15) 



assuming x(0) = y(0) = 0. Equations (TT4")) and (jT5]) show 
that cycloidal motion is possible in our system. This 
differs qualitatively from the Rashba s-o zitterbewegung 
which is always perpendicular to the initial v g . 

Figure 3 shows trajectories for three distinct v g = 
(hkoy / m*)y - all with ko y > 0. We find motion opposite 
to and along the y-axis (orbits I and III, resp.) and even 
a closed path (II). To understand this behavior we note 
that for nkoy << e — e e the linear-in-t terms in {ynit)} 
can be recast into hka y t/m*_ =>■ the injected wave moves 
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FIG. 3: Zitterbewegung due to the s-o coupling g for distinct 
ratios fi = e_/2e so . Note the peculiar trajectories with the 
forward injected electrons moving backward (I) and even in 
a closed path (II). This follows from the s-o induced change 
in the curvature of the bands which renormalizes the effective 
masses. Here we use \ so ko y = m/10, A^ 1 = m*r)/h. 



with the renormalized velocity v * = hko v /m*_ 
for orbit I: /j, < 1 ^ m*_ < 
II: / 



/' 



1 



1 



=>■ m_ 
ml > 



oo and v* = 
v* > 0. 



and 



Hence, 

and v* < 0; for orbit 
0; and for orbit III: 
Though remarkable, 



we stress that the orbits I and II occur for unusual pa- 
rameters (e.g., £f < e S o/10). However, these orbits do 
show that our s-o Hamiltonian has a physical mechanism 
allowing for cyclotronic motion without magnetic fields. 

Spin Hall conductivity <J z xy . The spin Hall effect is a 
convenient probe for s-o effects in wells [26[ ■ We have cal- 
culated a% y ( "clean limit" ) in the presence of an external 
magnetic field B by following the approach of Rashba 
271 ]. which allows us to properly account for both the 
intra- and inter-branch contributions in the Kubo for- 



mula 
find a 



xy 



Here we focus on the B — ► limit where we 
for sp > e (two subbands occupied) and 



for 



< e F 
2e so /e_ 



e 
< 



+ 



«2 + Ki/2 



(16) 



(upper subband empty), where 
= (ef — e+)/2e_, and K3 = 
\/ k1 /4 + K1K2 + 1/4. Note that a% y is non-zero and non- 
universal in this range, shows a discontinuity at Ef = e Qi 
and changes sign as a function of ep. Details of our calcu- 
lation of a * a nd a thorough discussion will be presented 
elsewhere |29| . Here we just note that measurements of 
a xy ( vs - £ f) m symmetric two-subband wells offer a pos- 
sibility to probe our new s-o interaction fjjj ]. Note that 
the Rashba and the linearized Dresselhaus spin Hall con- 
ductivities are identically zero in the DC limit 27 , 3(| . 

We have introduced an inter-subband-induced s-o in- 
teraction in quantum wells with two subbands. The cor- 
responding s-o coupling 7] (whose magnitude is similar 
to Rashba coupling) is non-zero even in symmetric wells. 
This new s-o interaction gives rise to a non-zero spin Hall 
conductivity, renormalizes the bulk mass by ~ 5% (mea- 
surable via cyclotron resonance [23| ) in double wells and 
can induce a cycloidal zitterbewegung. Weak antilocaliza- 
tion [2^, [3l[ should offer another possibility to measure 
77. 
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